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We develop a phenomenological thermodynamic theory of ferroelectric BaTi03 (BT) thin films 
cpitaxially grown on cubic substrates using the Landau-Devonshire eight-order potential. The con- 
structed "misfit-temperature" phase diagram is asymmetrical. We found that, overall view of the 
phase diagram depends on the values of compliances used in calculations and provide two qualita- 
tively different diagrams. A thermodynamic path for BT film grown onto a particular substrate can 
be found using a plot of the room-temperature tetragonal distortion (c — a) /a as a function of misfit 
strain. 
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I. INTRODUCTION 

Epitaxial ferroelectric BaTiC>3 (BT) films have been 
intensively investigated for their applications in high 
dielectric layers of capacitors, in dynamic random ac- 
cess memories (DRAM), microwave integrated circuits, 
and micro electro mechanical systems (MEMS) 0, 
Nanoscale ferroelectric heterostructures lead to new phe- 
nomena. The physical properties of ferroelectric thin 
films substantially differ from those of bulk ferroelectrics. 
Epitaxially grown thin films are usually highly con- 
strained because their fabrication is accompanied by 
several strain factors such as misfit strain due to lat- 
tice mismatch between the film and the substrate or 
the intermediate buffer layer, thermoelastic strain gen- 
erated by the difference between the thermal-expansion 
coefficients (TECs) of the film and the substrate, self- 
strain of the paraelectric-ferroelectric phase transforma- 
tion if the film is deposited above the Curie tempera- 
ture (T c ), and defects such as dislocations and vacan- 
cies, appeared during the deposition. The effects of 
internal stresses on the electrical and electromechan- 
ical properties have been investigated experimentally 
and theoretically for a number of ferroelectric systems 

a i a a a a a m ed, tn, m m ee m cli m tn. 

The experimental investigation of the above-mentioned 
factors is difficult because the physical properties of real 
thin films are a combined result of these factors. There- 
fore, theoretical analysis is important since it can provide 
fundamental insights into the behavior of thin films. 

A phenomenological thermodynamic theory of ferro- 
electric BT thin films epitaxially grown on cubic sub- 
strates was developed [5[ using the Landau-Devonshire 
six-order potential as a polynomial of the polarization 
components. For single-domain BaTiC>3 film, the "misfit- 
temperature" phase diagram was constructed. It was 



found that the two-dimensional (2D) clamping of the 
film, apart from a shift of the temperature of the fer- 
roelectric transition, changes the entire phase diagram, 
creating new phases that are not present in bulk crys- 
tals. For BT thin films epitaxially grown on cubic sub- 
strates Pertsev et al [5| have predicted five distinct stable 
phases. Four of them, namely p, c, aa and r have a com- 
mon multiphase point on the diagram at zero misfit strain 
and temperature equal to the T c of the bulk material. 

Later [f| , another phase diagram was proposed for BT 
film using another set of coefficients in the same six-order 
thermodynamic potential. In the vicinity of T c and small 
misfit strains, the phase diagram is quite similar to the 
previous one, except the low-temperature range where 
the 2D clamping stabilizes mainly the r-phase and not 
the ac one. Recent first principle calculations [TH, [H| 
qualitatively confirmed the latter finding and general 
view of the phase diagram. However, according to 
the multiphase point has a positive shift along the misfit 
strain axis. 

The misfit-temperature phase diagram constructed in 
Q, is based on the six-order thermodynamic poten- 
tial with the temperature-dependent coefficients [2(| Hl[ , 
which reproduces quite well the low-temperature ferro- 
electric phase transitions in BT single crystal. It is worth 
noting that the 2D clamping increases considerably the 
ferroelectric phase transition temperature in perovskite 
thin films. In this case, the coefficient at the sixth- 
order term in thermodynamic potential becomes nega- 
tive; therefore, the six-order expansion is not valid for 
high-temperature phase transitions in thin films. To 
overcome the problem, high-order terms should be in- 
cluded. In this paper we employed an eighth-order poly- 
nomial for the Landau-Devonshire potential of bulk BT 
single crystal [22| under mechanical stress-free boundary 
conditions, only the first coefficient at p 2 is temperature 



2 



dependent, while others are independent of temperature. 
We discuss temperature-misfit-strain phase diagrams for 
BT films on cubic substrates and their thermodynamic 
path determined using room-temperature tetragonal dis- 
tortion value. 



Starting from the Helmholtz potential F(j>, u) with me- 
chanical forces or stresses as external parameters, the fol- 
lowing expression is readily obtained for Gibbs potential 
<!>(p, a) = F(p, u) — au where u = u(a) can be found from 
the similar condition: 



II. A PHENOMENOLOGICAL 
THERMODYNAMIC POTENTIAL 

Following Pertsev et al. [j| we consider the thermo- 
dynamical phenomcnological potential for BT thin film 
epitaxially grown in a cubic paraelectric phase on a thick 
cubic substrate. Standard Gibbs potential $(p, a) is a 
function of polarization - p and stress - a . Consider- 
ing strain u as an external parameter (either geometrical 
parameters or stresses are defined) the Helmholtz poten- 
tial should be written as F(p,u) — <f>(p, a) + au, where 
a = cr(u) can be found from the following requirement 



da 



(1) 



dF(p, u) 
du 



(2) 



Equations (1) and (2) link the coefficients of the cor- 
responding potentials. In the absence of external forces 
and at p = 0, the strains u are governed by the thermal 
expansion. As follows from Eq. (1), the Gibbs poten- 
tial should contain the — chq{T — T a ) term, where ao is 
the TEC of the bulk material, T a is the temperature at 
which thermal deformations are absent. In the following 
we assume T a = 0. 

The Helmholtz phenomenological eighth-order poten- 
tial for cubic perovskite is given as [22TI24T]: 
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The last term in the F(p, u) expression (3) is respon- 
sible for the thermal expansion [23| of the cubic crystal. 
The strains Uj are given in the Voigt matrix notation. 
The film is under short-circuit boundary conditions. 



Here we consider (001) surface of a cubic substrate, 
which exhibits no symmetry change in the whole temper- 
ature range. The substrate defines the following strains 
at the film/substrate interface: tti = 1*2 = u s ,Uq = 0. 
They are governed by the symmetry of the substrate, 
the primary strain occurred at the film deposition, and 
TEC of the substrate. We assume the elastic fields to be 
homogeneous in the film volume and ignore any changes 
on the surface of the substrate. In the absence of external 
forces the remaining strains are determined by minimiza- 
tion of the Helmholtz potential (Eq. 3) with respect to 
the strains u,(« = 3, 4, 5), 
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Taking into account Eqs. 4 we can write the following 
thermodynamic potential of the thin film: 



G = a.\ (p\ 
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(5) 

The following relations between coefficients of poten- 
tial (5) and those in F(p,u) and <I>(p, a) are (note, coef- 
ficients used in terms Gq and G 8 are not changed under 
this transformation): 
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Equations (6) evidently contains linear thermal expan- 
sion of the film, so that misfit strain is u m = us — aoT 
Here us contains the primary film strain uo and ad- 
ditional strain due to thermal shrinkage of the substrate 
appeared on cooling from the deposition temperature. 
The primary deformation appears due to misfit of the 
lattice parameters of the film and the substrate at the 
deposition temperature To. Starting from the deposition 
film can be under stress if c (out-of plane) and a (in- 
plane) parameters of the film at To are different. Taking 
into account Eqs. (4) one can obtain the following ex- 
pression: 
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where c and a are the lattice parameters of the film at 
T . 

Shrinkage of the substrate on cooling from To down T 
is obvious, A = b \t — b \t = &o"s(T — To) , where as is 
TEC of the substrate and bo the lattice parameter of the 
substrate at Tq. Then, film strain in the plane parallel 
to the substrate is 



u s 



-an ao(l+aoTo+no)+A — ap 

bo , 



(8) 



= a To + Uo + ^as(T-To). 
Therefore, (u s — aoT) in Eqs. (6) is 

u rn =us — a Q T = u a - (a - —a s )(T - T ) (9) 
In the following calculations we apply this to Eqs. (6). 

III. PHASE DIAGRAM 

The phase diagram depends on the substrate used for 
film deposition because Eq. 9 contains lattice parameters 
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FIG. 1: Temperature-misfit strain phase diagram BT film for 
the eight-order thermodynamic potential with the coefficients 
taken from column I in Table f . Solid lines - first order phase 
transitions, dashed lines - second order phase transitions. The 
calculated coordinates of the multiphase points are T\ : (u m = 
0.01 xfO" 3 ,T = 388.2A"), T 2 : («m = 0.21 x I0~ 3 , T = 377K), 

u m = -0.39 x 

-3 



r 3 
10 



0.31 x 10" 3 ,T 



366K), T 4 

,T = 354RT), T 5 : (u m = -6.20 x 10^, T = 1UK) 
Thin lines 1 — 1 and 2 — 2 correspond to the thermodynamic 
path for BT film epitaxially grown with primary strain uo ~ 
on MgO and SrTiOa, respectively. 



and TEC of the substrate. In this case Eq. (9) deter- 
mines the line that is the thermodynamic path for partic- 
ular film/substrate heterostucture on the "temperature- 
misfit strain" (T — u m ) diagram. The slop of the ther- 
modynamic path on the diagram is driven by the sub- 
strate/film lattice parameters ratio and by the difference 
of their TECs, so that 

u m = u + (a -a s )T . 

This approach is valid if and only if coefficients at p 2 
are linear in T and u m . That is not true in the case of 
the sixth-order potential [13, Hl| . One can expect rather 



small changes on the phase diagram when (cvo 



bo 



as) 



is small. However, we found the following quantitative 
changes of the T — u m phase diagram: (i) The stability 
range of the ac phase is limited in the low-temperature 
range when the difference between TEC of the film and 
the substrate is small; (ii) a tricritical point appears on 
the line between paraelectric and ferroelectric c-phase. 

It is worth noting that values of the coefficients in Eq. 
(5), and therefore, overall view of the phase diagram de- 
pends on the values of elastic moduli or compliances used 
in the calculations. Elastic moduli in potential (3) should 
be taken at D = const, therefore, we used cfj|25| mea- 
sured at room temperature. The corresponding calcu- 
lated compliances are given in the second column of the 
Table 1. On the other hand, available compliances mea- 
sured ptjj , in the high-temperature paraelectric phase are 
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TABLE I: Coefficients of the thermodynamical potential (5) 
for BaTi03 , thin film calculated from the eighth-order poten- 
tial (3) with the coefficients taken from [23] a" = 4.124 x 
10 5 (T — 402), and two sets of compliances taken from [25j-I 
(calculations) and [2(|-II. 
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an 


2.56 
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5.97 


4.13 
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4.48 
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1.294 


1.294 


x 10 9 Jm 9 /C 6 
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-1.95 


-1.95 




ai23 


-2.5 


-2.5 
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—2.5 


—2.5 




aim 


3.863 


3.863 


x 10 10 Jm 13 /C 8 


"1112 


2.529 


2.529 




«1122 


1.637 


1.637 




«1123 


1.367 


1.367 




Qh 


0.10 


0.10 


m 4 /C 2 


Ql2 


-0.034 


-0.034 




Q44 


0.029 


0.029 




Sll 


10.79 


8.33 


xW- 12 m 2 /N 


S12 


-4.365 


-2.68 




S44 


7.94 


9.24 





"«i = 4.124 x 10 5 (T - 388 - 24.9 x 10 3 u m ) 
a a 3 = 4.124 x 10 5 (T - 388 - 28.3 x 10 3 « m ) 
s ai = 4.124 x 10 5 (T - 388 + 25.7 x 10 3 u m ) 
b a 3 = 4.124 X 10 5 (T - 388 + 29.2 X 10 3 u m ) in Jm/C 2 units 
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FIG. 2: Temperature-misfit strain phase diagram BT film 
for the eight-order thermodynamic potential with the co- 
efficients taken from column II in Table 1. Solid lines - 
first order phase transitions, dashed lines - second order 
phase transitions. The calculated coordinates of the mul- 
tiphase points are Ti : (u m = 0.0002 x 10~ 3 ,T = 388A"), 
T 2 : {u m = 0.20 x 10~ 3 ,T = 3737T). Thin lines 1-1 and 
2 — 2 correspond to the thermodynamic path for BT film 
epitaxially grown with primary strain uq — on MgO and 
SrTiOs, respectively. 



slightly different. We used coefficients for eighth-order 
potential from Ref. (22J. According to expression (6) 
two sets of elastic moduli produce two sets of coefficients 
of the thermodynamical potential for thin film listed in 
Table 1. Below we provide phase diagrams for two sets 
of coefficients. 

The temperature-misfit strain phase diagram of BT 
film for the eighth-order potential with the coefficients 
taken from column I in Table 1 is shown in Fig. 1. Among 
eight allowed by symmetry phases, only six are present 
on the diagram, 

p — p = (0, 0, 0)P4/mmm(i¥123) paraelectric, 

c p = (0, 0,p)P4mm(iV99), polarization normal to 
the film surface, 

a — p = (p,0,0)Pmm2(N25), polarization along the 
axis in plane of the film, 

aa — p = (p,p,0)Amm2(N38), polarization along the 
diagonal in plane of the film, 

r — p = (pi,pj,p2)Cm(N8), polarization with aa and 
c components [27| 

ac — p = (pi 1 0,p2)Pm(N6) 1 polarization with a and 
c components. 

The phase diagram for the eighth-order potential (Ta- 
ble 1, column I) contains one more phase - a with the 
polarization along one of the former cubic axis in plane 
of the film. This phase does not exist in the sixth-order 
Pertsev's model 6]. As shown in Fig.l, the a-phase is 
stable in the narrow range between paraelectric and or- 
thorhombic aa phase. Note, p a transition is second 
order, while a — aa is first order. Another important 
difference with respect to Q is the first order nature of 
the p — c transition with very small hysteresis. Also, in 
contrast to Pertsev's model [6[ the second order c r 
transition is possible only at very low temperature below 
T 5 point (u m = -6.20 x 10" 3 ,T = 1UK). Above T 5 
transition to the r-phase is possible through the inter- 
mediate ac-phase as a result of the second order c — > ac 
transition, and following the first order ac — > r one. As 
one can see in Fig. 1, the boundary between c and a 
phase in between T\ and T2 points, between c and aa in 
between Ti and T3, between ac and aa in between T3 
and T4 corresponds to the first order transitions. Note, 
transition between r and aa is always the first order. 
Finally, the T x point (u m = 0.01 x 10" 3 , T = 388.2iT) 
is slightly shifted with respect to zero misfit strain be- 
cause in contrast to Pertsev's model [6] p — c transition 
is of first order one. Such a positive shift along the mis- 
fit strain axis is in a qualitative agreement with the first 
principle calculations [l8l ]. 

The temperature-misfit strain phase diagram of the BT 
film for the eighth-order potential with the coefficients 
taken from column II in Table 1 is shown in Fig. 2. Two 
intermediate phases a and ac disappear and only four 
phases p, c, aa and r are present on the diagram. The 
p - c transition is of first order-type, but hysteresis is 
smaller. Therefore, shift of the T\ point with respect to 
zero deformation is now less pronounced. 

Thin lines 1 — 1 and 2 — 2 in Figs. 1 and 2 correspond 
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FIG. 3: Room-temperature tetragonal distortion (c — a) /a 
as a function of misfit strain for eight-order thermodynamic 
potential with the coefficients taken from column I in Table 
1. 

to the thermodynamic path ( Eq.(9) ) for BT film epitax- 
ially grown at Tq — 1023!^ with primary strain uq = 
on MgO and SrTi03, respectively. The primary strain 
Mo depends on the growing conditions and formation of 
misfit dislocations in the film volume. Real thermody- 
namic path is shifted with respect to 1 — 1 and 2 — 2 
lines on the actual uq value, which is determined from 
experiment. 

Actual Mo value can be found from the structural data. 
Figure 3 shows room-temperature tetragonal distortion 
(c — a) /a as a function of misfit strain u m calculated 
using the following expression: 



in Fig. 3, one can find room-temperature u m value that 
determines the corresponding point on the phase diagram 
(Figs. 1 or 2). The distance between this point and I — 1 
(2-2') line at T = 300K is equal to u - 



IV. CONCLUSIONS 

We develop phase diagrams for epitaxial BT films on 
cubic substrates as a function of the misfit strain based 
on the eight-order Landau-Devonshire phenomenological 
potential, which is valid for high-temperature states, fn 
the frame of this model the phase transition from the 
paraelectric to the ferroelectric c-phase is always of first- 
order-type. This result has important consequence, the 
phase diagram is no more symmetrical with respect to 
the zero misfit strain axis. 

Available from the literature, two sets of elastic moduli 
for bulk BT yields different phase diagrams with a dif- 
ferent number of stable phases. Therefore, phase transi- 
tion sequence in epitaxial thin films is governed not only 
by the misfit strain value u m . Other important factors 
are thermodynamic conditions which can change chem- 
ical composition of the film. In the case of ferroelec- 
tric perovskites, oxygen pressure during film deposition 
seems to be a very important factor, because even weak 
deviation of the film density markedly changes its elastic 
properties. 

The performed calculations allow to estimate possible 
phase transition temperatures for BT films on different 
substrates. We propose a simple way how to find a ther- 
modynamic path on the phase diagram for BT film grown 
on MgO or SrTi03 substrate using room-temperature 
structural data. 
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